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ABSTRACT 

By comparing the widths of absorption lines from O I, Si II and Fe II in the redshift Zabs = 
2.076 single-component damped Lya system in the spectrum of Q2206 — 199 we establish 
that these absorption lines arise in Warm Neutral Medium gas at ^ 12000 ± 3000 K. This 
is consistent with thermal equilibrium model estimates of ^ 8000 K for the Warm Neutral 
Medium in galaxies, but not with the presence of a significant cold component. It is also 
consistent with, but not required by, the absence of C II fine structure absorption in this system. 
Some possible implications concerning abundance estimates in narrow-line WNM absorbers 
are discussed. 
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1 INTRODUCTION 

The complex structure of the interstellar medium in the Galaxy 
may be considered as a mixture of four components: the neutral 
atomic hydrogen is found to be generally in cold neutral medium 
(CNM, < 500 K) and warm neutral medium (WN M, ~ 5000 
- 10000 K) gas iWolfe, Prochaska & GawTseill2003h . and ionized 
material may be warm (W IM, ~ lO" K) or hot (HIM, T ~ 10** K) 
dMcKee & Ostrikeilll977h . The WNM and CNM are expected to 
be in approximate pressure balance (e.g. Wolfire et al., 1995); HI 
at temperatures between the two is unstable and will rapidly evolve 
towards one of the two stable states. The CNM and WNM are usu- 
ally studied using HI 21cm absorption against background radio 
sources, and show that a significant fraction of the warm material 
may be in the unstable regime (see the discussion by Wolfire et al., 
1995, also Begum et al., 2010). 

Temperature estimates for Damped Lyman-alpha absorption 
systems (DLAs) reveal both CNM and WNM gas. The pres- 
ence of C II fine structure absorption is associated with the CNM 
I Wolfe. Prochaska & Gawiserl 20031; Wolfe. Gawiser & Prochaskal 



l2003h . and lHowk. Wolfe & Prochaskj j2005l) use C II and Si II fine 
structure lines to show that a syste m at redshift z = 4.2 contains 
CNM material. lEehner et alj ( |2008|) use an upper limit to the C II* 
column density in the z = 2.377 DLA to show that the neutral 
gas in this system must be warm. Modelling of a few weak Mg II 
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absorbers also suggests temperatures of ~ 10000 K lljones et al.l 
I2OIOI) . From 21cm spin temperature studies iKanekar et al.l ( 1200^ 
suggest that a significant fraction of the HI at redshifts z > 1.6 is 
warm. However, this conclusion is controversial as it requires as- 
sumptions about the temperature and density distribution of the gas 
across the 100 pc dimension s of the backgroun d radio sources that 
have not been demonstrated jWolfe et al.l201 ih . 

The few attempts to estimate temperatures from the line 
widths in DLAs have relied on the identification of narrow ab- 
sorption components, nor mally from CI, which show up clearly 
as distinct sharp features jlorgenson et al.l [20091 ; iTumlinson et al.l 
|2010| ;I Carswell et al.ll201ll) . The corresponding temperatures of a 
few hundred degrees Kelvin or less indicate that the absorbing gas 
is in a cold neutral medium (CNM). However, attempts to deter- 
mine temperatures by using the widths of singly ionized heavy 
element lines have not been very succes sful, either in the Ga laxy 
jNehme et al.l2008h or at higher redshifts jCarswell et alj201lh be- 
cause of uncertainties introduced by the blending of many velocity 
component s. One exception is for the interstellar medium near the 
Sun, where iRedfield & Linskvl l l2004t) find temperatures averaging 
6800K with a dispersion of ~ 1500 K. 

One way of avoiding velocity component confusion is to se- 
lect systems for which the velocity structure in the singly ionized 
species is simple. These are rare, and may not be typical, but it is of 
interest to see if the inferred temperatu re in such cases is ~ 8000 K, 
as would be expected for the WNM l l Wolfire et"al]|l995l) . or if a 
temperature of < 500 K typical of the CNM is more appropri- 
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ate. The other requirements are that at least two ion species with 
different atomic masses and transitions with a range of oscillator 
strengths be measured, and that at least some of these fall on differ- 
ent parts of the curve of growth (normally linear and logarithmic) 
so that th e Doppler parameters for the ions may be e stimated indi- 
vidually jjorgenson et al.ll2009lJCarswell et al jlOl ih . 

There are two studies where the inferred line widths have 
yielded temperature estimates in QSO absorbers dominated by 
single-component systems with strong H I, both as part of D/H 
abundance studies. lO'Meara et al.l ( 1200 lb use Keck HIRES spec- 
tra to derive a temperature T = 1.15 ± 0.02 x 10'' K, and a 
bulk motion component 6turb ~ 2.56 ± 0.12 km s^^ by com- 
paring the Doppler parameters for the HI, Nl and 01 lines in a 
system with HI column density logA'^(Hl) = 19.42 at redshift 
z = 2.536 towards the QSO HS0105 + 1619. The other is a DLA 
at redshift z = 2.076 in the spectrum of Q2206 - 199 which 
was studied byjP rochaska & Wolfe ( 1997) using Keck HIRES data. 
IPettini & Bowed {2001) used a combination of HST STIS and VLT 
UVES data to derive a similar, if less certain, temperature by com- 
paring the Doppler parameters for H I and singly ionized heavy el- 
ement lines. 

We revisit the second case here, since the temperature is based 
largely on low S/N STIS spectra for higher order Lyman lines. With 
the presence of several transitions of a range of ions of different 
masses, O I, Si II and Fe II there is the prospect of temperature esti- 
mation from the heavy element lines alone. In the next section we 
describe the spectroscopic data and the techniques used to establish 
a temperature T ~ 12000 ± 3000 K. Section [23l describes a series 
of tests undertaken to ensure that the result is reliable, and in partic- 
ular that the error range is realistic and that the absorbing medium 
must be warm, and some further temperature indicators are consid- 
ered in Section |2!4l Section [3] examines some implications of the 
temperature estimate, particularly for relative element abundance 
estimates in cases where the absorption lines are intrinsically nar- 
row, and the main conclusions are summarized in Section|4] 

2 THE REDSHIFT 2.076 ABSORPTION SYSTEM IN 
Q2206-199 

2.1 The data 

The absorption system at redshift z = 2.076 in Q2206 — 199 was 
fo und to have low heav y element abundances, with [Si/H] < —2.4, 
bv lRauch et"aD ( Il99(]l) using 30 km s ^ resolution spectra. Using 
higher resolution Keck HIRES spectra. IProchaska & Wolfel ( ll99'^ 
noted the very low Fell and AlII abundances in a system at red- 
shift 2 = 2.076214, and th at the system appeared to have only 
a single velocity component. IProchaska et alj 1 I2OOII) re-examined 
this system and further determined the abundance for Si II, further 
confirming it as a low abundance system with [Fe/H] = —2.6, and 
[Si/H] = -2.3. 

The spectra used here were obtain ed using UVES on the ESO 
VLT, for programmes 65.O-0 158(A) jPettini et alj|2002l) . 072.A- 
0346(A) jPetitiean et alj2006h and 074.A-0201(A). For each expo- 
sure ESO provides a seeing estimate from their differential image 
motion monitor. The standard setting DIC2 437 (3730 - 4990 A) 
& 860 (6650 - 10600 A) was used for two exposures for a total of 
7800s in 0.5 arcsec seeing as part of programme 65.O-0158(A). 
The other settings used in that programme were DIG I 390 (3260 - 
4450 A) & 564 (4580 - 6680 A), for 10800s with 0.5 arcsec seeing 
and 6300s with slit-limited (I arcsec) resolution. For programme 



072.A-0346(A) the seeing approximately matched the slit width of 
I arcsec, and the exposures totalled 9000s for DICl 346 (3030 

- 3880 A) & 580 (4760 - 6840 A) settings. Programme 074. A- 
0201(A) used a 0.9 arcsec slit in seeing ranging from 0.8 to 1.1 arc- 
sec with total exposure times of 19603s in settings DIC2 455 (3930 

- 5 160 A) & 850 (6580 - 10420 A). The data were reduced using 
the UVES pipeline and the UVES_popler packag^. The signal-to- 
noise ratio (S/N) in the coadded spectrum is 60 per 2.5 km s~^ 
pixel at 5000 A. 

The variety of slit widths and wavelength settings is likely to 
result in the spectral resolution being a function of wavelength. For 
individual exposures appropriate width Gaussians provide an ade- 
quate approximation to the instrument profile. For the slit-limited 
cases, which dominate the final sum here, we verified this assump- 
tion by establishing that single Gaussians gave g ood fits to n i eht sk y 
emission lines in the sky spectrum described bv' Hanuschi3(l2003l) . 
For the seeing-limited observations there may be extended wings 
beyond a core which is reasonably close to a Gaussian (see e.g. 
Diego 1985, King 1971). We have constructed an instrument pro- 
file as a sum of Gaussians with a common centre and with widths 
determined using the resolution - slit width product Rx s = 41400 
(s in arcsec) for the blue arm and Rx s = 38700 for the red arm of 
the spectrograph, weighting each by the relative efficiency (see the 
ESO UVES instrument characteristics webpag^ and the response 
curves linked to from there). For the 0.5 arcsec seeing-limited ex- 
posures we took s to be the seeing estimate. For the other expo- 
sures the resolution was close to being slit-limited, so s was set to 
the appropriate slit width. While for the seeing limited exposures 
the adopted profile will not give the full description, any flux out- 
side the core will be dominated by the slit-limited observations so 
should not be important. In any case, as we show below, the results 
here are insensitive to the line-spread-function adopted so a precise 
one is not needed. 

The error estimates provided by the data extraction package 
were found to be close to those inferred by examining the root- 
mean-square (RMS) deviations in continuum regions. However, in 
lines which have zero residual intensity the RMS measures indi- 
cated that the error estimates were too low by a factor of 2. We ver- 
ified that this is generally true for UVES spectra extracted in this 
way by examining data for other objects. To correct for the error un- 
derestimate when the data is near zero we multiplied the error term 

by (1.0 + 15.0 X min (l.O, max (0.0, (1.0 - di/c^)f) ) ^ , where 
di, Ci are the data and continuum values in the zth data pixel. While 
it is not clear what the true correction should be, this function was 
chosen because it changes the error little when the data is not close 
to zero and doubles it near zero intensity. Provided that the correc- 
tion is sharply peaked near the data zero, the precise form of the 
correction makes little difference to the results here. 



2.2 Temperature estimation 

The full optical wavelength coverage of these spectra allows us 
to measure several Fell and Sill lines with a range of oscillator 
strengths / and wavelengths A, and Ol 1039 and 1302 with /A 
values differing by 0.8 dex. So, at least in principle, we may esti- 
mate the Doppler parameters for these species independently. Also, 
with a range of /A values the profile fits depend most strongly on 
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Figure 1. The Voigt profile fits and data for the C II, O I, Si II & Fe II lines used on a velocity scale relative to redshift z = 2.0762286. The data are shown as 
black histograms and the overall fits in blue. Where they differ from the overall fitted profiles, the fits to the transitions of interest here are indicated are shown 
in red, and any interlopers shown in green. All but Sill 1808 are shown against a unit continuum, and the continua for different fines have been separated by 
an integer offset. The Si II 1808 line is weak, and has been scaled up by a factor of five so it can be seen more clearly. The dip (in red) at +26 km s^^ relative 
to the Sill 1260 fine is Fell 1260 at the same redshift. 



the curve of growth, and not on an estimate of the instrument pro- 
file jCarswell et al.ll2oTll) . We assume that all ions arise predomi- 
nantly in the same region where a single temperature is appropriate. 
Since they arise in a DLA these assumptions should be close to re- 
ality. Within the region we assume that the Doppler parameter for 
a given ion is 6 = i/^turb + 2fcT/m, where &turb is the turbulent 
(bulk) component, k Boltzmann's constant, m the ion mass and 
T the temperature. We used the VPFIT prograr rjf] with oscillator 
strengths from the compilation bv lMortonI ( |2003|) . For the Doppler 
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parameters fcturb and T were used as the independent variables for 
at least two ions of different mass at the same redshift simulta- 
neously, with the constraints that both variables are non-negative. 
The error estimates from the program then apply to feturb and T, 
not to the Doppler parameters for the individual ions. The temper- 
ature component of the Doppler width varies as l/y'm, so even 
for O I and Fe II the thermal widths differ by less than a factor of 
two and for other available ions the ratio is even smaller. If there 
is a significant turbulent component then the differences in the line 
widths could be quite small. Under these circumstances, even quite 
well-determined Doppler parameters for individual ions can yield 
large error estimates for both 6turb and T. The Doppler parameter 
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Table 1. Q2206 - 199 z = 2.076 system profile fit parameters. 
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Column densities N have units cm ^ , and Doppler parameters h are in km s ^ . 

'??' denotes an unidentified line, taken to have the same rest wavelength and oscillator strength as Lya. 



error estimates shown in parentheses in Table[T]are given as a guide 
only, and were obtained by assuming that the 6- value for each ion 
was determined independently of the others. 

The redshifts of CII 1036 and 1334 were tied to be the same 
as for the other ions to provide a (small) increase in the precision 
of the z-estimate. However, since the C II /A values differ by only 
0.14 dex they were not used for constraining the temperature. The 
Nl lines at 1200 A are too weak to be useful. 

In the z = 2.076 absorption system Sill 1260, 1304 and 1808 
are reasonably clear of contaminating absorption lines from sys- 
tems at other redshifts, but Sill 1526, which is clear of the Lya 
forest, is unusable because it is badly blended with a strong com- 
plex Fell 1608 at redshift z — 1.9205. Allowance has to be made 
for blended Lya forest lines for all usable Sill lines except 1808. 
This is straightforward generally, since the Lya lines are signifi- 
cantly broader than the heavy element ones. Some of the longer 
wavelength Fe II lines are affected by atmospheric absorption and 
emission, but Fell 1608, 2374, 2382 and 2600 are clear of serious 
contamination. At shorter wavelengths Fell 1260 is in the wing 
of Sill 1260, and Fell 1143 and 1144 provide useful constraints. 
For Ol 1302 and 1039 are also useful, though the 1039 line is 
blended with a broad Lya. The results of fitting these with VP- 
FIT are given in Table [T] along with the blends fitted at the same 
time. The Doppler parameters for 01, Sill and Fell were deter- 
mined using common values of bturb and T, which were treated as 
independent variables. The eiTor estimates for the Doppler param- 
eters given in brackets are derived assuming that they are indepen- 
dent of those for the other ions. The line designated by '??', which 
is assigned the same wavelength and oscillator strength as Lya, is 
a marginal narrow component which is probably an unidentified 
heavy element. 

There are two data uncertainties which should be allowed for. 
The continuum level cannot be determined precisely, especially in 
the Lya forest, and inaccurate continuum levels can cause inconsis- 
tencies when several lines of the same ion are fitted simultaneously. 



Table 2. Other z = 2.076 system ions 
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z 
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0.01 
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0.06 
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2.0753686 


0.0000121 


7.54 


1.88 


12.09 


0.08 


CIV 


2.0761144 


0.0000058 


5.11 


1.49 


12.67 


0.20 


CIV 


2.0762127 


0.0000286 


17.32 


2.09 


13.57 


0.04 


CIV 


2.0764230 


0.0000136 


9.57 


1.87 


13.00 


0.22 


CIV 


2.0766579 


0.0000409 


12.42 


5.11 


12.24 


0.17 


SilV 


2.0761218 


0.0000020 


7.75 


0.33 


12.64 


0.02 


SilV 


2.0764017 


0.0000065 


12.91 


1.05 


12.55 


0.03 



Note: H I, C II, Al II and Al III redshifts and 6- values are obtained from the 
parameters given in Table [T] 



The other is the true zero level for the data, which can affect strong 
lines particularly, can differ from the data extraction zero by 1 or 
2%. The VPFIT package allows both of these to be used as further 
free parameters in the overall fit, and we have done so where nec- 
essary. We do not give their values in Table [T] since they depend 
on how well the continuum and zero levels have been estimated 
originally. 

From Table[T]we see that the temperature estimate for the sys- 
tem is r = 1.22 ± 0.32 X 10" K, so consistent with it being WNM 
but not CNM. The turbulent component of 5.3 km s~^ is larger than 
the thermal broadening for any of the ions at the inferred temper- 
ature. For 01 the thermal width is 6 — 3.5 km s^^, for Sill 2.7 
km s~^ and for Fell 1.9 km s^^. However, these thermal widths 
are enough to give measurable differences in the Doppler parame- 
ters for the ions involved. 
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The CII column density given in Table [T] is the best fit ob- 
tained by leaving its Doppler b to be determined independently of 
the other ions. If we impose the turbulent component and temper- 
ature derived above, then for CII 6 = 6.75km s~^ and log A*" = 
14.38 ± 0.03 (see Tabled. 

For the fitted profiles given in TableEthe reduced = 1.103 
(X^ = 893.39 for 810 de grees of freedom). The formal probability 
that this fit describes the data is 0.022. One might improve the fit by 
adding more blended lines, but the the reduction in is to some 
extent offset by the loss in the number of degrees of freedom - see 
iRauch et al. for a discussion on this point. We have not quite 

reached that stage: the temperature estimate does depend on hav- 
ing identified all the narrow (< 10 km s~^) blends with key lines. 
For example, adding an extra narrow component to Si II 1304 gives 
a system temperature of 11800K. However, the Sill 1260 profile 
suggests that this narrow component is not Sill, and we have not 
found any other lines in a system which would plausibly identify 
such a component. 

For this syst em the HI column d ensity is logA^(Hl) = 
20.43 ± 0.06 jPettini et al.l 1 19941 : IProchaska et all 1200 lb . 
IPettini et all ( 1 1994 show a fit to the Lya profile. From the UVES 
data we obtain a very similar result (log A'^(Hl) = 20.44 ± 0.05, 
with the error dominated by systematic effects). However, it is 
difficult to be sure of the continuum over the damped hya pr ofile 
in the echelle spectrum so we adopt the IPettini et alj i 19941) and 
IProchaska et al.l ( l200lh value. 

That the singly ionized species do arise in a neutral region, so 
with the temperature established above, in the WNM, is strongly 
indicated by the fact that it is the only component which is likely to 
give rise to the damped Lya line. Further, the corresponding Al III 
lines are weak, and from Table |2]the AlIII/AlII ratio is —0.7 dex, 
so the singly ionized species dominates. There is also no sign of a 
significant component of Si IV at the Si II redshift, though it would 
fall between the two Si IV components given in the table. There is a 
C IV component at a redshift which, within the rather large errors, 
is consistent with the low ionization redshift. However the Doppler 
parameter is high so it is unlikely to be closely associated with the 
region of interest. 

If the ionization correction is negligible for the low ioniza- 
tion co mponent, then the a bundances relative to the solar photo- 
sphere l|As2iundetai]200i) are [Fe/H] = -2.59±0.06, [Si/H] = 
-2.29 ± 0.06 and [O/H] = -2.28 ± 0.08. Using the turbulent- 
plus-thermal value for the Doppler parameters, [C/H] = —2.64 ± 
0.06, so [C/0] = —0.36. These abun dances and c o lumn densi- 
ties are consistent with those given by IPettini et~aL I ( l2008l) . who 
assumed a common Doppler parameter of b = 6.5km for all 
ions. The Fe II and Si II values are also similar to those reported fo r 
this system by IProchaska & Wolf el ( ll997l) : |Prochaska et al.l ( 1200 ll) . 

2.3 Is the result reliable? 

The temperature estimate here is 12200K, but the VP- 
FIT error estimate of 3200K is large, so at the la level 
9 X 10^ < T < 1.54 X lO^K. When the errors are roughly com- 
parable with the best estimate value one might be concerned at the 
reliability of the error estimate, since e.g. negative temperatures are 
unphysical. This is not quite the case here, but to further check we 
have generated artificial spectra using the best-fit parameters given 
in Table[T] convolved these spectra with the appropriate instrument 
profiles, added noise to mimic that in the original data, and de- 
termined the fit parameters from these spectra using the same fit 
regions and allowance for continuum and zero level adjustments 




5000 10"* LSxIO"* 2x10* 

temperature (K) 



Figure 2. Voigt profile fit values as a function of temperature from fits 
to 1, Si II & Fe II lines. The dashed lines show the Ax^ level for the la 
range, and the dotted lines for 2(t, for two degrees of freedom. The . 
value is 893.39 (811 d.f.). 



as we did for the original data. From 2000 trials we find that the 
mean fitted temperature is 12270 K, with a standard deviation of 
3350 ±55 K. The expected error in the mean temperature estimated 
from the 2000 trials is then 3350/V1999 = 75 K. This mean tem- 
perature is therefore not significantly different from the input value, 
while the standard deviation is a little higher than the mean VPFIT 
estimate of 3200 K. 

Another way of checking the results is to derive the best fit 
X^ values for a range of assumed temperatures T, and then use 
the values of Ax^{T) = X^{T) — Xmin to estimate the error 
ranges . The method is described bv iLampton. Margon & BowveJ 
( Il976l) . In the application here we have at least two closely corre- 
lated (or anti-correlated) variables, the temperature and the turbu- 
lent Doppler parameter, so for a given probability P we look for 
the temperature range for which Ax^(T) is less than the value of 
X^ corresponding to P for two degrees of freedom. Under this as- 
sumption, for a degree of confidence of 68.3% (Ict for a normal 
distribution) Ax^iT) < 2.3, and for a < 2a deviation for a two- 
sided threshold) Ax^iT) < 6.2. The results are shown in Fig.|2] 
from which it is evident that the best fit value for the temperature is 
in good agreement with that obtained from VPFIT, but the la error 
estimates from this method, ±3900 K, are a little pessimistic, with 
those from this procedure being a little over ~ 15% larger than 
those from the simulations. 

These analyses preclude the possibility that a significant 
amount of the absorption arises in cold gas. If we take 1000 K 
as a somewhat generous upper limit for CNM gas, then using the 
Doppler parameters are almost but not quite the same for the three 
ions, we find Ax^ = 18.86. The formal probability of this occur- 
ring by chance is 1.4 x 10~^. 

Another way of establishing the probability that the material 
is WNM rather than CNM is to generate synthetic spectra with the 
best-fit parameters for an assumed temperature of 1000 K, which 
we take as the upper limit for the CNM. These were log A'^(0 1) = 
15.23 (cm-2) with 6(Ol) = 5.99km s'\ logiV(Sill) = 13.69, 
6(SiIl) = 5.95km s~^ and log iV (Fell) = 13.36, fo(Sill) = 
5.93km s^^. As in the simulations described above, we convolved 
these synthetic spectra with the appropriate instrument profiles, 
added noise to mimic that in the original data, and determined 
the fit parameters from these spectra using the same fit regions 
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and allowance for continuum and zero level adjustments as we did 
for the original data. We also did this for 2000 trials, and find in 
only one case was the fitted temperature above 10000 K, with none 
above UOOOK. Since our temperature estimate was 12200 K, we 
conclude that the probability that we are observing CNM gas is 
< 0.0005. 

Since we have chosen lines which sample different parts of 
the curves of growth we do not expect that our results will depend 
strongly on the chosen instrument profile. However, since the ab- 
sorption lines are resolved there could be some dependence on its 
assumed shape and width. We have undertaken some tests and find 
that this is not an important consideration. These were: 

(i) As a first check we note that for lines in the rest wavelength 
range 1140 - 1350 A and for Sill 1808 the seeing limited (~ 0.5 
arcsec) and slit limited (~ 1 arcsec) exposure times were compara- 
ble, in the range 1.2:1 to 0.7:1. For Fell 1608, the Fell lines above 
2300 A, and O I 1039 the ratio is ~ 0.5: 1 or less. Since the inferred 
Doppler parameter for Fe II is the lowest and most of the lines used 
in its determination have the broadest instrument profiles, it sug- 
gests that the raw instrument profile is not propagating directly into 
the inferred line widths. 

(ii) We have further checked that the temperature estimate is 
largely independent of the assumed instrument profile by perform- 
ing the Voigt profile fits with different weights for the 0.5: 1 arcsec 
components. We also used single Gaussian point-spread functions 
for the whole spectrum, with widths ranging from 4.0 - 8.0 km s~^ 
FWHM. While the details changed, and at the ends of the Gaussian 
ranges the fits were not really acceptable, the temperature estimate 
proved robust with 11000 < T < 13000 K. The error estimates 
were similar in all cases, at < 3600 K. 

While the narrow lines in the Lya forest are easily separable 
from the broad Lyman lines, the effective continuum for these nar- 
row lines is affected by the Lya absorption. The fitting procedure 
does take account of the additional uncertainties given the fit model, 
but, since the O I 1039 line is strongly blended with Lya, as a rough 
consistency check we have refitted using only Sill and Fell to es- 
timate the temperature. Not surprisingly, since the total ion mass 
range is lower by almost a factor of two, the error estimates in- 
creased significantly. We obtained 6turb ~ 5.2 ± 0.3km s~^and 
T = 1.7 ± 0.7 X lO" K. So the best fit is at a higher temperature, 
but without the O I the significance level for the material to be in 
the WNM rather than the CNM falls to 2.3cr. 

Ultimately, the most compelling evidence for a thermal com- 
ponent could come from curve-of-growth considerations, since 
then there are no uncertainties involving the instrument profile. If 
there is a thermal broadening component, then different mass ions 
which have observed transitions with a good range of /A values 
will lie on different curves of growth to fit the observed equivalent 
widths. The compatibility of VPFIT and the cur ve of growth analy- 
ses has been illustrated bv lCarswell et alj ( l201 ih . However they are 
not likely to be completely equivalent. Changing from a description 
which for a set of lines assigns a single redshift, two parameters to 
describe the Doppler widths, and individual column densities for 
each ion to one where separate single parameters (the equivalent 
width) are derived for each line and then combining them as curves 
of growth will not necessarily yield exactly the same results. Also, 
blends with other lines are involved as in the Lya forest, the indi- 
vidual equivalent widths are more difficult to determine. 

We have estimated the equivalent widths for the lines used by 
adding the contributions from all pixels wholly within ±20 km s~^ 
of the line centres, with error estimates obtained by summing the 




log(NfX) 

Figure 3. Curves of growth using the fitted Doppler parameters given in 
Table[T]for OI (purple), Sill (blue) and Fell (red), showing the estimated 
equivalent widths for the lines used. Upper limits ai'e shown for blended 
lines, along with best estimates (with errors) after the blends were removed. 
The horizontal error bars are based on the column density uncertainties from 
Voigt profile fitting. The lower panel gives detail near the logarithmic parts 
of the curves of growth. See text for details. The units for NfX are cm^ ^ . 



pixel contribution errors in quadrature. In this case possible blend- 
ing is particularly important, since only lines with rest wavelengths 
> 1410 A are clear of the Lya forest, and so both O I lines are pos- 
sibly blended and for Si II only the line at 1 808 A is clear. Fe II has 
several unblended features with a range of /A values, but the tran- 
sitions at 1260 A and shortward are blended. So for the curves of 
growth based on the Doppler parameters given in Table [T] shown 
in Fig. [5] several of the lines appear as upper limits where they 
are blended. For these we have also estimated equivalent widths 
over the ±20 km s~^ range after removing extraneous blends us- 
ing profiles generated from the values given in that table. These 
are shown as points near the curves underneath the limit markers. 
These curves of growth do not cover all the lines used in the Voigt 
profile fits. Fell 1260 is close to Sill 1260, and so was omitted here. 
Also Fell 1144 is close to a strong absorption feature, so while its 
upper limit is shown (at log NfX = 7.33) its corrected equivalent 
width was subject to large errors so it was omitted. The point at 
logNfX = 7.32 near the curve corresponds to Fell 2374, which 
is free from contaminating blends. 

While Fig. [3] illustrates that the profile fitting results are con- 
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sistent with a curve of growth analysis with different curves for dif- 
ferent mass ions, the curves of growth alone do not require different 
curves for the various ions. They allow both fully turbulent and the 
mixed turbulent-plus-thermal broadening solution given above. We 
suspect this broader tolerance range comes from the broader error 
ranges allowed by this independent line approach. This is improved 
a little by changing the pixel weighting by fitting independent Voigt 
profiles to each component and then estimating equivalent widths 
from the line parameters, but that did not reduce the errors by very 
much. To obtain a significant result we require stronger constraints 
e.g. by having a common Doppler parameter and column density 
for all lines of a given ion, but then we lose the independent er- 
rors so the curve-of-growth picture is less useful. Curves of growth 
based directly on the fitted ions were, not surprisingly, compatible 
only with some WNM thermal broadening component. 

Fig. [3] also illustrates another point which should be high- 
lighted. Where the lines are narrow, a reliable Doppler parameter 
determination for any ion requires a number of transitions with 
a range of /A values which cover different parts of the curve of 
growth. At the least not all the transitions should be only on the 
logarithmic or the linear parts of the curve since in either case the 
Doppler parameter is poorly constrained. Also, if ions with single 
transitions (or transitions with similar fX values) are available they 
should not generally be included in the joint turbulent/thermal fit. 
Since the column density for each is then a free parameter, they pro- 
vide minimal additional information on the intrinsic line widths. 
Where the lines are marginally resolved the results may even be 
misleading if such lines are included. We illustrate this by remov- 
ing Ol, for which two lines with fX values separated by 0.82 dex 
were used, from the temperature analysis (Table [T]l and try using 
C II in its stead. Si II and Fe II together do not constrain the temper- 
ature very well, so C II could, in principle, provide some improve- 
ment. However, since the two CII lines at 1036 and 1334Ahave 
similar fX values the result is critically dependent on the assumed 
line spread function (LSF). The nominal result for our best estimate 
of the instrumental profile is T = 24000 ± 3500 K. If we choose 
an instrumental FWHM for both CII lines of 5.5 km s~^ then the 
best fit temperature is T = 2.44 ± 0.36 x lO" K, while a FWHM 
of 8.0km s^^ gives T = 9.6 ± 4.8 x 10^ K (in both cases the error 
estimates are those given by VPFIT). So, unlike the case where 01 
was used, the result obtained with C II is dominated by uncertain- 
ties in the instrument profile. 

The results above are based on a single component model with 
a single temperature and Gaussian turbulent velocity distribution. 
Given the apparent differences in line widths between the various 
ions of different atomic masses it is hard to see how any mul- 
ticomponent model could give significantly different temperature 
estimates. We did attempt two- and three-component models con- 
strained to have only turbulent broadening in each velocity com- 
ponent, and each was a significantly poorer fit to the data than the 
model with some thermal broadening. 

We also tried a two-component turbulent-plus-thermal model 
and, not surprisingly, one component still had significant thermal 
broadening while another was completely turbulent. Because two 
close velocity components provided the best fit in this model there 
is considerable degeneracy between them, and so the error esti- 
mates in individual quantities are so large as to make them almost 
meaningless. 

A possibility is that some of the O I absorption might arise in 
a fully neutral region which is distinct from the region giving rise 
to the other lines, and then any small velocity shift between the two 
might be interpreted as a broadening of the O I lines. However, this 



is unlikely since there is no detectable redshift difference between 
the 01 and other species, and nor is there any sign of CI which 
would also arise in such a neutral component. 



2.4 Other temperature indications 

IPettini & BowenI ( 120011) describe an analysis to determine the 
D/H ratio in this system using a spectrum obtained with STIS 
on the HST covering the high order Lyman lines from HI 926 
down beyond the Lyman limit. This spectrum covers the wave- 
length range 2775 - 2860 A with a continuum S/N~ 6 per 9.5 
km s^^ pixel. Since the Lya line alone gives the HI column den- 
sity from the damping wings, and the high order Lyman lines will 
fall on the logarithmic part of the curve of growth, it should be 
possible to estimate the HI Doppler parameter assuming that the 
absorption arises in a single component. With hydrogen the mass 
range is much greater than is available with heavy elements alone, 
so despite the poor S/N and resolution it should further constrain 
the temperature estimate. IPettini & BowenI ( 1200 ih used VPFIT on 
these data and estimated a temperature of ~ 11000 K assuming a 
thermal component of 5 km s^^. This is in extremely good agree- 
ment with the value we have derived, but we suspect part of this 
agreement is fortuitous. Not only is the S/N poor in the STIS spec- 
trum so blended Lya at lower redshifts are difficult to identify, but 
also the continuum is very uncertain, the zero level not well de- 
termined (it is about 4% too high in the Lyman limit region), and 
the line profiles not well sampled. From joint fits Including both 
the STIS and UVES data, with the assumed HI column density 
logiV(Hl) = 20.43, we find that 5.2 < fcturb < 5.4km s"^ and 
9.6 X 10^ < T < 13.6 X 10^ K. Much of the range comes from so- 
lutions involving different zero level and blending assumptions. We 
have not explored these very far, so the parameter ranges should be 
taken as indicative rather than firm, but they suggest that we have 
gained some precision but not a large amount. The deuterium col- 
umn density logA''(Dl) lies in the range 15.2 - 15.9 (cm~^), so 
the possible error ra nge for D/H is la r ger (b y a factor of almost 
3) than suggested b y IPettini & BowenI ( 1200 1). This reinforces the 
IPettini et al.l ( l2008h suggestion that the earlier error estimates had 
been too small. 

There is another approach which can yield temperature lim- 
its to DLA absorbe rs, so we have explored its applicability here. 
iLehner et alj ^20081) have used the CII*/HI column density ratio 
in a component of a redshif t z ~ 2.4 DLA and the models of 
IWolfe. Prochaska & Gawise3 ( l2003h to show that the neutral gas 
in that system must be warm. The reason in that case is that the 
CII* 1335 is so weak relative to HI Lya that CII 158/im cannot 
be a dominant coolant so, unless the star formation rate is low, the 
absorption must be in WNM gas. Since we have shown that the 
gas in the case studied here is WNM the C II*/H I is also likely to 
be low. The upper limit to the C II* column density in this system 
is logiV(ClI*)< 12.75 (2cr, TablelU and with logiV(Hl)= 20.44 
the C II 158 fim cooling rate is log Ic < —27.2 erg s^^ (H atom)~^. 
This limit is too high by about 0.4 dex to require a WNM, but is 
consistent with it. 
3 IMPLICATIONS 

Some previous temperature estimates at high redshift have been 
based on CI and provided uppe r limits of a few hundred Kelvin 
djorgenson et al.l2009l ; ICarswell et al. 2011 ), so ind i cative of CNM 
gas. WNM material was found by ,0'Meara et alj ( 1200 ih using a 
long atomic mass baseline from HI to effectively Sill. Here we 
have shown that it is possible to use a range of ions from O I to 
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Fell to differentiate between cold and warm neutral media. A re- 
quirement is that the absorption lines used for each ion have a range 
of /A values sufficient to constrain their curves of growth, and if 
this condition is satisfied for all ions used then the results will not 
be strongly dependent on the instrumental point-spread-function. 
So in favourable cases it should be possible to obtain temperature 
estimates without the large mass range provided by hydrogen and 
deuterium with the heavy elements. However, the lines which are 
used should not be strongly blended, and the turbulent component 
should not predominate. The h ~ 5 km s~^ found in the case de- 
scribed here is probably close to the upper limit for current echelle 
data. Suitable cases may still be quite rare, since most of those in 
the literature exhibit complex line profiles. A preliminary examina- 
tion of over 20 UVES spectra of DLAs indicates that for > 90% of 
them the components are too badly blended to consider further. 

Much of the emphasis on high redshift absorption systems 
has been to examine the redshift evoluti on of the heavy e lemen t 
abundances, and most of these, e.g. by IProchaska et alj l l2003h . 
have relied mainly on the apparent optical depth method (AODM) 
dSavage & S embach 1991; Sembach & Savage 1992). This method 
generally uses transitions which are near the linear part of the curve 
of growth, and so the results are insensitive to the Doppler param- 
eter or any velocity structure. For DLA gas, abundant ions such 
as Sin and Fell have transitions with a wide range of oscillator 
strengths so their column densities can be estimated reliably. For 
less abundant species such as Znll, CrII and Nill the lines are al- 
most invariably sufficiently weak that the method gives satisfactory 
results. In principle the O I column density could be constrained by 
its lower oscillator-strength transitions, but these fall in the Lya 
forest, and so they are subject to greater uncertainties. 

In the case of the Zabs = 2.076 DLA in the Q2206- 199 spec- 
trum the heavy element column densities are not sensitive to the 
line broadening model adopted since bulk motion is a significant 
contributor to the line widths for all ions. However, this may not 
always be the case. For an example of how the assumed broaden- 
ing function may affect the inferred column densities we consider 
the cas e of CII at z = 2.340 towards J0035 - 0918 l lCooke et al.l 
l20Ilah . They found a Doppler parameter b = 2.36k m s~^ which 
they us ed for all ions. For the data described by ICooke et al.l 
( 1201 lij) only Ol and Si II have transitions with a range of /A val- 
ues, and so these two species were used. From these, the mate - 
rial may be purely turbulent, as assumed by ICookeetalJ feOIlah . 
purely thermal, or anything in between with a best fit temperature 
of T ~ 4000 K. If instead of a purely turbulent model we con- 
sider the other extreme of a purely thermal model (i.e. turbulent 
component assumed to be zero) applied to those data, the temper- 
ature is T = 7.66 ± 0.57 x 10^ K. Using this thermal model the 
column densities of the heavier ions are close to the values given 
by Cooke et al., but now log 7V(C II) = 14.47 ± 0.06 (instead of 
15.47 ± 0.15 for a CNM model) and logiV(Ol) = 14.73 ± 0.06 
(CNM: 14.96 ± 0.08). The very large difference in the CII value 
arises because the C II lines are saturated and on the logarithmic 
parts of the relevant curves of growth. Then [C/Fe]= 0.51 ± 0.10 
rather than 1.53, which is still marginally outside the nominal range 
—0.1 < [C/Fe] < 0. 4 found in other low abundance and very 
high redshift syst ems jPenprase et al.l I2OI0I ; ICooke et alj 1201 Ibl : 
iBecker et al]|2012h . 

While the example above is an extreme one, it does reinforce 
the need for care in inferring column densities, as o pposed to lower 
limits as in e.g. iDessauges-Zavadskv et alj l l200ll) . where only a 
few saturated lines are available. 



4 CONCLUSIONS 

Comparison of the line widths for a range of ions of differ- 
ent masses, 01, Sill and Fell, in the Zabs = 2.076 system in 
Q2206-199 yields a temperature estimate of 12000±3000 K. This 
is consistent with expectations of ~ 8000 K for the Warm Neutral 
Medium in galaxies, but not with the material being predominantly 
in a Cold Neutral Medium. In this case there is a significant bulk 
motion common to all ions, which is normally modelled as turbu- 
lence. Independent of the nature of the bulk motions, it is difficult 
to avoid the conclusion that the intrinsic line widths increase as the 
atomic mass decreases, and that there is a significant component 
for which temperatures are ~ 10* K. 

Previous temperature estimates have relied on the comparison 
of H I and D I Do ppler parameters with those from heavy elements 
dO'Meara et al.l2001 ; Pettini & Bowen 2001) . We have shown here 
that if the ions used cover an adequate mass range, e.g. O I to Fe II, 
and for each ion the transitions cover a wide enough range of os- 
cillator strengths that different parts of their curves of growth are 
adequately sampled, then useful temperature estimates may be ob- 
tained. If these conditions are satisfied, then the temperature will 
not generally depend strongly on the assumed line-spread-function 
(LSF) for the spectrograph. However, including other ions where 
the conditions are not satisfied can lead to the result which depends 
strongly on the assumed LSF. 

A significant fraction of DLA gas may exhibit temperatures of 
~ 10* K. This is predicted f or gas with weak C II* absorption b y 
two-phase models for DLAs jWolfe. Proch aska & Gawiser 2003), 
and may be indicated by the H 1 21 cm results ( Kanekar et al. 2009]). 
In such cases the inferred column densities for important species 
such as C II, where the two observable interstellar lines have similar 
oscillator strengths (or more strictly, /A values) can depend criti- 
cally on the assumed broadening mechanism for narrow-lined sys- 
te ms. For example, the very high carbon overabundance reported 
bv lCookeet"aLl 12011 ah could be reduced by a factor of ~ 10, from 
~ 35 to ~ 3, if the lines are purely thermally broadened. 

Finally, we restate something which is well-known but some- 
times overlooked: the column densities for a given ion are uncertain 
and depend on the assumed Doppler parameter only when all the 
transitions for that ion are saturated and when thermal broadening 
component could be significant. In cases where the Doppler param- 
eters give implausibly high temperatures t hen e .g. single curve of 
growth analyses such as that by iMortod ( 1 19781) should give reli- 
able results for the ion column densities. The column densities in 
the z — 2.076228 abs orber reported here are also consistent with 
those reported earlier l lProchaska & Wolfelll997l : I Prochaska et al.l 
[200 1; Pettini et al. 20(51) because there is significant bulk motion 
so a single Doppler parameter is a reasonably good approxima- 
tion. However if the Doppler parameters for the heavier metals 
are < 3 km s~^ then thermal broadening may be important, and 
should be allowed for in estimating column densites of abundant 
species such as C II and O I. It is unfortunate that C II, a species of 
current interest in low metallicity cases, provides only two ground- 
state transitions for interstellar medium absorption studies. These 
are usually saturated and, if that is the case, any derived column 
densities are likely to be uncertain. 
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